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SUPERSONIC  FLUTTER  OF  A  THERMALLY  STRESSED  FLAT  PANEL 


WITH  UNIFORM  EDGE  LOADS* 

By  Harry  G.  Schaeffer  and  Walter  L.  Heard^  Jr. 
Langley  Research  Center 


SUMMARY 


The  effect  of  a  nonunlfonn^  self-equilibrating  stress  state  associated 
with  a  parabolic  temperature  distribution  is  included  in  the  flutter  analysis 
of  a  rectangular^  sirr^ly  supported  panel  subjected  to  airflow  over  one  surface 
and  uniform  edge  loads.  Llneap?^zed^  static^  two-dimensional  aerodynamics  is 
used  in  conjunction  with  thiii-plate  theory  to  formulate  the  problem,  y^j^umerical 
results  for  a  square  panel  show  that  nonuniform  self -equilibrating  stresses  are 
as  important  as  uniform  stresses  in  affecting  flutter  behavior.'^  ^  ^ 


IMTRODUCTIOH 


Dynamic  instability  of  thin  platelike  structures  in  a  gas  flow  is  a  prob¬ 
lem  which  has  received  considerable  attention  in  recent  years.  As  a  result  of 
the  high  velocities  attainable  by  modern  aircraft^  severe  aerodynamic  loading 
of  exterior  skin  panels  frequently  occurs  and  can  lead  to  panel  flutter.  (See, 
for  example,  ref.  1.)  In  the  treatment  of  panel-flutter  problems,  theoretical 
investigators  have  taken  into  account  such  effects  as  panel  configuration,  flow 
orientation,  boundary  conditions,  and  uniform  midplane  stress.  However,  effects 
of  thermally  induced,  self-equilibrating,  midplane  stress  have  not  been 
examined . 

It  has  already  been  shown  in  reference  2  that  a  system  of  uniform  stresses 
can  greatly  reduce  the  flutter  speed  of  an  unbuckled  panel.  The  purpose  of  the 
present  paper  is  to  investigate  the  effects  of  a  nonuniform,  self-equilibrating, 
stress  state  on  the  flutter  behavior  of  a  flat,  uniformly  stressed  panel.  Since 
a  nonlinear  temperature  distribution  could  easily  result  from  aerodynamic 
effects,  it  is  assumed  for  the  present  analysis  that  the  nonuniform  stress  state 
is  thermally  induced  by  a  parabolic  temperature  distribution  which  is  constant 
through  the  thickness.  However,  the  results  are  valid  for  any  similar  system 
of  self -equilibrating  midplane  stresses. 


*Part  of  the  information  presented  herein  was  offered  as  a  thesis  by  the 
second  author  in  partial  fulfillment  of  the  requirements  for  the  degree  of 
Master  of  Science  in  Engineering  Mechanics,  Virginia  Polytechnic  Institute, 
Blacksburg,  Virglna,  December  1964. 


For  the  purpose  of  this  analysis  the  panel  is  considered  to  be  simply 
supported  against  lateral  deflections,  but  unrestrained  to  inplane  displace¬ 
ments  at  the  boundaries.  In  order  to  simplify  the  analysis  the  aerodynamic 
forces  are  assumed  to  be  those  resulting  from  Ackeret  theory.  The  analysis 
proceeds  by  using  the  Galerkin  technique  in  conjunction  with  an  assumed  stress 
function  to  determine  an  approximate  midplane  stress  distribution  corresponding 
to  the  ass\amed  temperature  variation  and  boundary  conditions.  This  thermal 
stress  distribution  is  then  substituted  into  the  linear  partial  differential 
equation  governing  lateral  deflections  of  thin  plates  in  combined  bending  and 
compression.  The  Galerkin  procedure  is  used  again  to  obtain  approximate  solu¬ 
tions  for  the  flutter  behavior,  and  numerical  results  are  presented  for  a 
square  panel  subjected  to  various  combinations  of  uniform  loading  applied  at 
the  boundaries-. 
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SYMBOLS 

amplitude  coefficients  (see  eq.  (ll)) 
panel  length 
panel  width 
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panel  stiffness. 
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E  Young ' s  modulus 

h  panel  thickness 

n^  s  integers 
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s^(^2)rs^(^3)rs  integrals  defined  by  equations  (ijb)  to  (ijd) 
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(ff  ^  (f  I 


'rs 
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aerodynamic  loading  coefficient 


Lrs 


Z  aerodynamic  pressure  load  per  unit  area 

M  Mach  number 

Mrs  coefficient  defined  by  equation  (A5) 

Nx,Ny,Nxy  midplane  stress  resultants 

KxQ^Ny-Q  midplane  uniform  stress  resultants 

Nx^,Ny.^,Nxy^  midplane  thermal  stress  resultants 


^rs 

q. 

Qrs 
R 

Rxo  =  — i 


coefficient  defined  by  equation  (A6) 
dynamic  pressure 

coefficient  defined  by  equation  (A7) 
total  number  of  terms  in  flow  direction 
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S 
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T 


total  number  of  terms  in  cross-flow  direction 
time 

temperature 

maximum  amplitude  of  temperature  distribution 


lateral  deflection  of  panel 
Cartesian  coordinates  of  panel 
coefficient  of  thermal  expansion 
mass  density  of  panel 
when  i  =  j 
when  i  /  J 


A 

^cr 

cp 


2qa5 

dynamic-pressure  parameter^  — — r= 

dVm^  -  1 


critical  value  of  A 
Poisson *s  ratio 
Airy  stress  function 

p 

oEha 

thermal  stress  parameter^  - - - 


ijr  =  _ 
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n  generalized  eigenvalue 


m  frequency 


Notation 

=  (  ),xx  (  ),yy 
=  (  )^xxxx  ),xxyy  ^  (  ),yyyy 

When  subscripts  follow  a  comma^  they  indicate  partial  differentiation 
of  the  principal  symbol  with  respect  to  the  subscripts. 


STATEMEM'  OF  PROBLEM  AXD  ANALYSIS 


A  flat^  rectangular  panel  of  length  a^  width  b^  and  uniform  thickness  h^ 
as  shown  in  figure  1^  is  considered  in  this  analysis.  The  panel  is  simply  sup¬ 
ported  on  all  edges^  exposed  to  supersonic  airflow  at  Mach  number  M  over  one 
surface j  and  subjected  to  a  nonlinear  temperature  distribution  T  =  T(x_,y)  in 
the  middle  plane.  The  panel  edges  are  considered  to  be  unrestrained  to  thermal 
expansion;  however^  uniform  loading,  Nxq  and  which  is  independent  of 

the  temperature  distribution  may  be  applied  at  the  boundary. 


Basic  Equations 

In  this  analysis  thin-plate  theory  including  the  effects  of  inplane  loads 
is  used.  The  basic  equilibrium  equations  for  this  case  are 


DV^ 
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(1) 


(2) 
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It  is  assumed  that  Nx  and  Ny  are  positive  in  compression  and  that  all  mate¬ 
rial  properties  are  to  he  evaluated  at  the  average  panel  temperature. 


An  approximate  expression  for  the  aerodynamic  pressure  load  I  is  obtained 
by  use  of  linearized,  static,  two-dimensional  aerodynamics  (Ackeret  theory). 

This  approximation  not  only  greatly  simplifies  the  analysis,  but  as  shown  in 
reference  5,  yields  essentially  the  same  flutter  results  as  exact  aerodynamics 
at  M  ^  1.5  for  uniformly  stressed  panels,  except  when  the  stress  level  causes 
two  modes  of  oscillation  to  have  the  same.  Or  nearly  the  same,  frequency.  Thus, 
if  care  is  taken  to  avoid  these  stress  levels,  Ackeret  theory  should  also  be 
applicable  to  the  nonuniform  stress  problem  of  the  present  analysis,  so  that 
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-2l 
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(3) 


The  midplane  force  intensity  terms  in  equations  (l)  and  (2)  are  written  as 
the  sum  of  those  induced  by  applied^  uniform^  normal  forces  at  the  boundary 
(subscript  o)  and  those  resulting  from  the  nonlinear  temperature  distribution 
(subscript  T)  as  follows: 


Ny  =  Ny^  +  N 


Xrp 


%  =  %o  %T 


Nxy  -  Nxy^ 


The  thermally  induced  stresses  are  determined  in  terms  of  a  stress  func¬ 
tion  cp  =  cp(x,y)  defined  by  the  following  relations: 


%T  -  'f. 


XX 


(5) 


Nxy,p  - 


This  stress  function  Identically  satisfies  the  equilibrium  conditions  given  by 
equations  (2).  For  compatibility  of  inplane  strains,  the  stress  function  must 
satisfy  the  following  partial  differential  equation  (ref.  k): 

V^^cp  =  aEhV%  (6) 

where  a  is  the  coefficient  of  thermal  expansion  at  the  average  plate  tempera¬ 
ture  and  the  temperature  T  is  a  function  of  x  and  y.  Note  that  the  sign 
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of  the  right-hand  side  is  opposite  to  that  normally  encountered^  because  direct 
stresses  are  positive  in  compression. 


The  panel  boundary  conditions  have  been  chosen  in  such  a  "way  that  they  can 
be  satisfied  by  using  rather  simple  functions.  Thus^  the  resulting  flutter 
analysis  is  greatly  simplified,  but  still  predicts  trends  which  one  may  expect 
to  be  valid  for  the  actual  panel. 


The  condition  that  the  panel  be  free  from  thermally  induced  normal  and 
shear  stresses  on  the  boundaries  req,uires  that  the  stress  function  satisfy  the 
following  boundary  conditions: 


<?(0,y)  =  q>(a,y)  =  cp(x,0)  =  cp(x,b)  =  0 


(7) 


The  condition  that  the  panel  boundaries  be  simply  supported  requires  that 
the  displacement  function  satisfy  the  following  boundary  conditions: 


w(O^y^t)  =  w(a_,y,t)  =  w(x,0,t)  =  w(x,b,t)  =  0 

(O^y^'t)  =  “  W^yy(x^0,t)  =  W^yy(x,b,t)  =  0 


(8) 
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The  analysis  which  follows  consists  of  two  parts.  First,  the  coefficient  of  an 
assumed  stress  function  cp  is  determined  approximately  by  using  the  Galerkin 
technique;  then,  the  stability  determinant  is  determined  from  the  deflection 
equation,  again  by  using  the  Galerkin  technique. 


Solution  of  Stress  Function  Equation 

The  temperature  distribution,  considered  in  this  problem,  varies  paraboli- 
cally  in  the  flow  and  cross-flow  directions  and  may  be  represented  by 

T  .  lto,(|)(l  -  |)(X)(l  -  f) 

where  is  the  maxlmiim  amplitude  of  the  temperature  distribution. 

In  order  to  simplify  the  subsequent  flutter  analysis,  it  is  assiomed  that 
an  approximate  stress  function  is  of  the  following  form: 

,  =  -  if  (9) 

This  approximation  of  the  stress  function  satisfies  all  the  boundary  conditions 
given  by  equations  (7)^  but  does  not  satisfy  equation  (6).  The  best  fit  is  made 
by  applying  the  Galerkin  technique  to  equation  (6)  to  evaluate  the  constant  C. 
The  result  is 


(10) 


The  accuracy  of  equation  (9)  in  approximating  the  true  stress  state  associated 
with  the  temperature  distribution  will  be  assessed  subsequently. 


Solution  of  Differential  Equation 

It  is  assumed  that  the  solution  of  equation  (l)  can  be  represented  as 
follows: 

.(x,y,t)  .  Ha  J  £  sin  ^  sin  2=  a^”"  (11) 

m=l  n=l 

vhere  the  frequency  cjd  is^  in  general^  complex.  Each  term  of  the  assumed 
solution  satisfies  the  given  boundary  conditions  on  w. 

After  substituting  equations  (j)^  (4),  (5)^  (lO)^  (H)  into  equa¬ 

tion  (l)  and  applying  the  Galerkin  procedure^  the  following  set  of  equations  is 
obtained  for  the  amplitude  coefficients  8^^* 


where 
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The  integrals  (li)^g^  (l2)^s^  (^3)rs  arise  due  to  the  nonuniform 

stress  distribution  and  the  integral  which  arises  due  to  aerodynamic 

loading  are  evaluated  in  appendix  A,  By  using  the  results  of  appendix  equa¬ 
tion  (12)  becomes 
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The  coefficients  Mrs^  '^rs  defined  by  equations  (Al), 

(A5),  (a6),  and  (AJ),  respectively. 

The  parameter  \|f  is  analogous  to  the  nondimens  tonal,  midplane,  uniform 
stress  parameters  Rxq  and  Ry^  and  may  be  thought  of  as  a  thermal  stress 

parameter.  The  range  of  the  indices  r  and  s  is  R  and  2S-1,  respectively, 
where  R  is  the  total  number  of  sine  terms  in  the  flow  direction  and  S  is 
the  total  number  of  sine  terms  in  the  cross-flow  direction.  Only  the  sine 
terms  corresponding  to  an  odd  number  of  half  waves  in  the  cross-flow  direction 
need  to  be  considered,  since,  for  the  symmetric  temperature  distribution  used, 
even  and  odd  cross-flow  terms  uncouple  and  the  odd  terms  lead  to  the  lowest 
(critical)  flutter  speed.  Thus,  equations  (l4)  represent  RS  (R-times-S) 
linear,  homogeneous,  algebraic  equations  for  the  RS  unknown  amplitude 
coefficients  a^s-  In  order  that  nontrivial  solutions  of  the  system  of  equa¬ 
tions  (l4)  exist,  it  is  necessary  that  the  determinant  of  coefficients  vanish. 
Thus,  the  stability  criterion  may  be  written  in  the  following  form; 

det^Aij  -  =  0  (15) 

vhere  is  the  eigenvalue. 

Since  the  problem  is  one  of  determining  the  stability  of  a  given  form  of 
solution^  it  is  most  natural  and  advantageous  to  associate  the  eigenvalue  with 
the  frequency  parameter  Then^  for  the  nonconservative  problem,  the  system 

is  dynamically  unstable  when  the  eigenvalue  0,  becomes  complex.  This  implies, 
according  to  equation  (ll),  that  the  system  diverges  in  an  oscillating  fashion. 
The  problem  now  is  to  determine  the  relationship  that  must  exist  between  the 
parameters  a/b,  ^Xq^  ^Yo^  ^  order  that  the  system  be  stable 

in  the  sense  described. 

For  the  case  where  R  =  2  and  S  =  1,  that  is. 


TCX 

ail  sin  —  +  api 


equation  (15)  simplifies  to  the  point  where  it  is  possible  to  determine  an 
analytical  expression  for  values  of  A  corresponding  to  harmonic  motion  of  the 
panel.  This  result  is 
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(16) 


The  critical  value  of  A  can  be  found  by  maximizing  A,  expression  (l6), 
with  respect  to  k2.  Solving  for  k2  and  substituting  the  result  into  eq.ua- 
tion  (l6)  gives 


‘cr 


9/ 

16 


(17) 


where  Acr  is  the  critical  value  of  the  dynamic-pressure  parameter  which  leads 
to  flutter. 


For  higher  ordered  modal  solutions,  the  algebraic  difficulties  in  obtaining 
expressions  similar  to  equations  (l6)  and  (l?)  are  formidable,  and  recourse  is 
made  to  an  iterative  procedure  using  a  digital  con5)uter.  For  the  computer  solu¬ 
tion  the  eigenvalues  are  calculated,  and  the  lowest  value  of  A  for  which  two 

of  the  eigenvalues  coalesce  is  sought.  The  approximate  modal  solutions  con¬ 
sidered  are  for  the  cases  where  R  =  4-,  S  =  1;  R  =  6,  S  =  1}  R  =  6,  S  =  2; 

and  R  =  6,  S  =  J.  The  coefficients  associated  with  the  latter  case  have  been 

calculated  from  equations  (l4),  and  the  general  matrix  equation  is  presented  in 
appendix  B.  The  coefficients  associated  with  the  other  cases  can  be  determined 
by  deleting  appropriate  rows  and  columns. 


DISCUSSION  OF  RESULTS 


Thermal  Stress  Distribution 

An  approximate  stress  distribution  has  been  assumed  in  order  to  simplify 
the  subsequent  dynamic  analysis.  Since  the  purpose  of  this  paper  is  to  deter¬ 
mine  the  effect  of  the  stress  state  associated  with  a  nonlinear  temperature  dis¬ 
tribution  on  the  flutter  characteristics,  it  is  important  to  consider  the  accu¬ 
racy  with  which  the  stresses  obtained  from  the  assumed  stress  function  approxi¬ 
mate  the  true  stress  state.  The  true  stress  state  associated  with  the  parabolic 
tenqperature  distribution  has  been  determined  by  the  finite-difference  method  of 
reference  5* 

Comparisons  between  the  one-term  approximation  and  a  finite-difference 
solution  are  presented  in  figures  2  and  J.  Figure  2  shows  the  variation  of 
Nx^  with  X  for  specified  values  of  y,  and  figure  3  shows  the  variation  of 

Nxyy  with  x  for  specified  values  of  y.  The  normal  stress  determined  by  the 
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one-term  approximation  and  that  determined  by  the  finite-difference  solution 
agree  quite  -well  (fig.  2).  In  the  case  of  the  shear  stress  distribution 
(fig.  5),  the  agreement  in  magnitude  is  not  as  close  as  that  for  the  direct 
stresses;  however,  it  appears  that  the  agreement  is  sufficiently  close  so  that 
the  overall  results  can  be  interpreted  with  confidence. 


Effect  of  Temperature  Distribution  on  Flutter  Boundary 

The  effect  of  the  parabolic  temperature  distribution  on  the  flutter 
behavior  of  a  square  panel  free  of  uniform  midplane  stress  is  shown  in  figure  4. 
The  results  are  presented  for  a  6-by-3  (R=6;  S=3)  term  solution.  The  inter¬ 
section  of  the  flutter  boundary  with  the  A-axis  is  the  value  of  Acr  associated 

with  the  unstressed,  unheated  panel.  This  value  decreases  as  i|f  increases 
(ATi  increases ).  until  the  flutter  boundary  becomes  tangent  to  the  curve  labeled 

"thermal  buckling  loop . " 

The  three  regions  shown  in  figure  4  are  characterized  by  the  value  of  the 
frequency  parameter  squared  (k^) .  In  the  region  labeled  "panel  flat,  no  flut¬ 
ter,"  k^  is  real  and  positive;  hence,  co  is  real.  In  the  region  labeled 
"flutter, "  k^  is  complex;  thus,  co  is  complex  and  at  least  one  root  will  lead 
to  oscillating,  divergent  panel  motion.  In  the  region  labeled  "panel  buckled, 
no  flutter,"  k^  is  negative;  thus,  co  is  pure  imaginary  and  the  panel  is 
statically  unstable. 

These  three  regions  are  separated  by  two  boundaries.  The  first  is  the 
buckling  loop  which  is  the  locus  of  points  for  which  k^  =  0.  The  second  is 
the  flutter  boundary  which  is  the  locus  of  points  at  which  two  frequencies 
coalesce.  The  flutter  boundary  is  terminated  at  its  point  of  tangency  with  the 
static  buckling  loop,  since  this  point  represents  the  limit  of  linear  plate 
theory.  This  transition  point  is  of  considerable  importance,  because  experience 
has  shown  that  it  represents  the  lowest  value  of  Acp  associated  with  this 

panel  configuration. 

Figure  4  shows  that  the  stress  distribution  associated  with  a  parabolic 
temperature  distribution  can  cause  a  6l-percent  reduction  in  the  value  of  Acr 
associated  with  an  unheated,  unstressed  panel.  In  order  to  give  an  indication 
of  the  magnitude  of  ATp  which  causes  an  effect  of  this  order,  ATp  is  com¬ 
puted  for  a  representative  square  aluminum  panel  with  ^  =  300.  The  results 

show  that  a  temperature  difference  of  only  27°  F  between  the  center  and  the 
edges  of  the  panel  causes  the  6l-percent  reduction  in  Acr* 

The  effect  of  the  parabolic  temperature  distribution  on  the  flutter  bound¬ 
ary  associated  with  unifomn  compressive  loads  in  the  flow  direction  is  shown  in 
figure  5  for  the  6-by-3  term  analysis.  The  curve  labeled  "\|/  =  0"  is  essentially 
the  exact  solution  established  in  reference  2  and  shows  how  the  critical  value 
of  the  dynamic -pres sure  parameter  varies  with  applied  uniform  loading  in  the 
flow  direction.  Increasing  values  of  t  lower  the  flutter  boundary. 
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The  regions  of  stability  and  instability  are  not  shown  in  figure  5  but  are 
similar  to  those  shown  in  figure  k;  however^  in  figure  5  the  dynamic -pres sure 
parameter  is  plotted  against  Rxq  instead,  of  t-  For  every  curve  represented 

by  some  constant  value  of  there  is  an  associated  buckling  loop.  The  dashed 
line  represents  the  locus  of  points  where  the  flutter  boundary  becomes  tangent 
to  its  corresponding  buckling  loop. 

The  decrease  in  due  to  a  parabolic  temperature  distribution  is  of  the 

same  order  of  magnitude  as  the  decrease  in  Acr  a  uniform  compressive 

load  in  the  x-direction.  The  comparable  reductions  in  Acr  "to  the  two  dif¬ 

ferent  parameters  indicate  that  the  effect  of  the  nonuniform  temperature  dis¬ 
tribution  is  as  important  as  that  of  uniform  compressive  stress.  Since  the 
nonuniform  stress  distribution  could  have  resulted  from  causes  other  than  a 
temperature  distribution^  the  previous  statement  can  be  generalized  to  say  that 
any  system  of  nonuniform  self-eq.uilibrating  stresses  may  have  a  significant 
effect  on  panel  flutter  behavior. 


Reference  2  shows  that  uniform  compressive  loads  applied  perpendicular  to 
the  direction  of  airflow  have  virtually  no  effect  on  the  flutter  of  a  uniformly 
stressed  panel.  It  is  of  interest,  therefore,  to  investigate  the  effect  of 
Ry^  for  the  present  case.  In  figure  6  the  relation  between  Acr  Fxq 

presented  for  various  values  of  the  stress  ratio  differences  in 

the  flutter  boundaries  due  to  Ry^  are  not  significant  and  cannot  be  plotted 

on  the  scale  presented.  However,  the  transition  point  between  the  flat-panel 
and  the  buckled-panel  flutter  boundaries  is  significantly  affected  as  indicated 
by  the  loci  of  termination  points  of  the  flat-panel  flutter  boundaries.  Thus, 
for  the  given  range  of  parameters  and  the  assiamed  nonuniform  stress  distribu¬ 
tion,  the  results  of  figure  5  may  be  used  to  predict  the  flutter  behavior  of  a 
square  panel  with  good  accuracy  for  any  combination  of  Rxq^  ^  which 


does  not  cause  the  panel  to  buckle.  Only  the  terminal  points  need  to  be 
located.  The  loci  of  these  terminal  points  for  several  values  of  the  stress 
ratio  are  shown  in  figure  6.  The  curves  apply  for  the  stress-ratio  range 


0  < 


lo 

xo 


<  1 


only,  since  large  positive  values  and  large  negative  values  of  the 


ratio  may  yield  different  results  because  of  changes  in  buckling  mode  shapes  in 

Ry 

the  cross-flow  direction.  When  —2.  =  the  locus  of  termination  points  of 

-^o 

the  flat-panel  flutter  boundary  is  essentially  a  horizontal  line.  Thus  for 
this  stress  ratio,  the  theory  predicts  that  the  transition  point  is  essentially 
independent  of  the  relative  magnitudes  of  applied  uniform  stress  and  thermal 
stress. 


Since  the  results  presented  in  this  paper  were  obtained  by  applying  the 
Galerkin  technique,  there  may  be  some  question  as  to  convergence  toward  the 
exact  solution  of  the  problem.  In  order  to  investigate  the  question  of  con¬ 
vergence,  results  for  various  approximate  solutions  are  presented  in  figure  7* 
The  flutter  boundaries  of  Aq^  as  a  function  of  Rx^  for  values  of 
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\|f  of  10,  20,  50,  and  40  are  presented  in  figure  7  for  the  2-by-l,  4-by-l, 

6-by-l,  6-by-2,  and  6-by-5  term  analyses.  The  solid  curves  represent  results 
of  the  6-by-2  and  6-by-5  term  analyses  whereas  the  other  curves  represent  the 
results  of  the  lower  order  solutions. 

An  examination  of  figure  7  shows  that  the  flutter  boundary  is  significantly 
altered  when  the  number  of  terms  in  the  flow  direction  is  Increased  from  two  to 
four.  However,  a  further  Increase  to  six  terms  exhibits  very  little  effect; 
thus,  the  solution  is  assumed  to  be  converged  for  flow- wise  terms.  The  require¬ 
ment  of  relatively  few  flow- wise  terms  for  convergence  is  not  unexpected,  since 
for  this  range  of  parameters  a  four-term  solution  gave  a  close  approximation  to 
the  exact  solution  in  reference  2  where  ^  is  zero. 

Inclusion  of  two  cross-flow  terms  in  the  analysis  instead  of  only  one 
exhibits  only  a  very  slight  lowering  of  the  flutter  boundary  as  may  be  seen  by 
comparison  of  the  6-by-l  and  6-by-2  solutions.  (For  the  scale  shown,  the 
6-by-l  and  4-by-l  resxilts  are  coincident  except  in  fig.  7(a)  where  the  6-by-l 
solution  gives  the  higher  values  of  Acr-)  Inclusion  of  three  cross-flow  terms 
gives  essentially  the  same  results  as  the  6-by-2  solution  and  indicates  that 
the  6-by-2  solution  is  converged  for  the  range  of  t  considered. 


CONCLUDING  REMARKS 


(llonuniform  stresses  associated  with  a  parabolic  temperature  distribution 
are  shown  to  be  as  Important  in  affecting  panel  flutter  behavior  as  uniform 
loading  applied  in  the  direction  of  airflow.  It  is  also  shown  that  although 
uniform  loading  applied  perpendicular  to  the  direction  of  airflow  has  virtually 
no  effect  on  the  flutter  boundaries  for  values  of  the  stress  ratio  Ryo/^Xo 

from  0  to  1,  the  transition  point  between  the  flat-panel  flutter  boundary  and 
the  buckled-panel  flutter  boundaiy  is  significantly  affected.  As  this  stress 
ratio  approaches  the  value  of  1,  the  transition  point  becomes  essentially  inde¬ 
pendent  of  relative  magnitudes  of  applied  uniform  loading  and  thermal  stress. 

Because  of  the  effects  of  the  nonuniform  stresses  shown  herein,  it  appears 
that  thermally  induced  nonuniform  stresses  should  be  considered  in  the  correla¬ 
tion  of  experimental  flutter  results  with  theory.  In  fact,  consideration  should 
be  given  to  any  system  of  nonuniform,  self-equilibrating  stresses,  even  those 
resulting  from  causes  other  than  a  nonlinear  temperature  distribution.  I  ^ 
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APPENDIX  A 


LOAD  COEFFICIEOTS 


Evaluation  of  the  integral  resulting  from  the  aerodynamic  loading 
(eq.  (13a))  gives  the  following  expression: 
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where  the  notation  m  /  r  indicates  that  the  terms  in  the  sum  associated  with 
m  =  r  are  to  he  omitted. 

The  integrals  (^1)^3^  (^2)^3^  (^3)rs  equations  (l3t>)^ 

(I3c)^  and  (l5d)^  respectively,  are  evaluated  by  making  use  of  the  specific 
definition  of  the  stress  function  given  by  equation  (9)-  The  result  is 
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As  in  equation  (Al)^  the  notation  m  ^  r,  n  ^  s  indicates  that  the  terms  in 
the  sums  associated  vith  m  =  r  and  n  =  s  are  to  be  omitted. 
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MATRIX  EQUATION  FOR  THE  UNKNOWN  AMPLITUDE  COEFFICIENTS 
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ONE-TERM  APPROXIMATION 
FINITE-DIFFERENCE  SOLUTION 


FLUTTER 


Figure  4.-  Thermal  flutter  boundary  for  square  panel  with  no  externally  applied  forces. 
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